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ABSTRACT 


This  report  describes  the  results  of  research  to  correlate  the 
fracture  resistance  with  microstructural  features  of  Salem  limestone  and 
Berea  sandstone.  Tests  were  conducted  on  wedge-loaded  double -cantilever -beam 
specimens  containing  machined  slots  as  crack  starters.  The  fracture  resistance 
of  Salem  limestone  measured  in  terms  of  R,  the  energy  dissipated  per  unit  area 
of  projected  surface,  was  found  to  increase  in  the  initial  stage  of  crack  ex¬ 
tension  but  finally  reached  a  constant  value  which  was  strongly  dependent  on 
crack  orientation  with  respect  to  the  bedding  plane.  For  this  limestone  R 
ranges  from  about  50  joules/m  to  230  joules/m  .  The  sandstone  which  is  a 
softer  and  more  compliant  rock  exhibited  similar  fracture  behavior  while  the  R 
values  ranged  from  465  joules/m  to  1580  joules/m  .  In  addition,  tests  in 
liquid  nitrogen  which  were  aimed  at  eliminating  plastic  deformation  in  the 
rocks  during  fracture  showed  little  difference  in  R  for  the  limestone  but  a 
substantial  reduction  in  R  for  the  sandstone  which,  in  the  latter  case,  may  be 
caused  by  differential  thermal  expansion  between  the  quartz  grains  and  the  cal- 
cite  cement.  Acoustic  emissions  were  detected  in  both  rocks  at  very  early 
stages  in  the  tests  indicating  the  occurrence  of  microcracking  near  the  initial 
slot  tip  at  low  loads.  The  mode  of  fracture  and  the  fracture  path  in  both 
materials  was  identified  by  fractography  over  a  broad  range  of  magnifications. 
The  evidence  gained  from  this  work  strongly  points  to  the  existence  of  an 
extensive  array  of  microcracks  produced  in  a  region  surrounding  the  main  crack 
tip  and  which  advances  with  it.  The  energy  dissipated  in  fracturing  of  rock  is 
associated  with  the  creation  of  the  large  amount  of  surface  area  contained  in 
this  microcrack  array.  This  picture  provides  a  self-consistent  explanation  for 
the  puzzling  discrepancy  between  the  measured  tensile  strength  of  rocks  and  the 
strengths  predicted  from  measured  R  values. 
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ON  FRACTURE  PROPAGATION  IN  ROCK 

by 

R.  G.  Hoagland 

G.  T.  Hahn,  Principal  Investigator 
A.  R.  Rosenfield 
R.  Simon 
G.  D.  Nicholson 

INTRODUCTION 

There  are  a  number  of  problems  which  arise  in  developing  relations 
between  the  fracture  behavior  of  rocks  and  their  mic restructure.  Among  these 
are: 

(®)  Variability:  The  composition  and  microstructure  of 
even  the  most  common  rock  types  are  usually  highly 
variable.  Hence,  mechanical  property  data  obtained 
from  small  laboratory  specimens  cannot  be  safely 
applied  in  the  field  if  one  disregards  the  possible 
scatter  of  the  data  due  to  microstructural  variations. 

(b)  High  Fracture  Energies:  It  is  often  observed  that  the 
mechanical  properties  of  rocks  are  considerably  differ¬ 
ent  from  those  of  the  individual  components  of  the  rock. 

This  is  particularly  true  in  regard  to  the  response  to 
a  tensile  stress.  The  energy  dissipated  in  propagating 
a  crack  in  rocits  is  typically  much  greater  than  is  con¬ 
sumed  for  the  single  crystals  of  its  constituents. 


Strengths  of  Smooth  Bars:  When  the  energy  dissipated  in 
propagating  a  large  crack  is  used  to  predict  tensile 
strengths  of  smooth  specimens  by  means  of  the  Griffith 
theory,  the  predictions  are  very  poor,  the  predicted 
strengths  greatly  overestimating  the  measured  strength. 
One  is,  therefore,  faced  with  the  puzzling  result  that 
rocks  are  apparently  much  tougher  (capable  of  dissi¬ 
pating  more  energy)  when  they  contain  a  large  flaw  than 
when  they  contain  a  small  one. 

(d)  Size  Effect:  There  is  often  a  significant  size  effect 
since  the  behavior  of  small-scale  laboratory  specimens 
or  models  may  not  compare  well  with  the  behavior  of 
large  in  situ  rock  masses.  Therefore,  to  ensure  that 
a  laboratory  test  is  to  provide  reliable  information 
-  is  necessary  to  recognize  the  important  microstuctural 

features  and  their  dimensions  and  to  scale  the  test  speci¬ 
fic  * 

mens  accordingly.  ' 

Comparison  with  Other  Brittle  Materials;  Rocks  can  be 
considered  an  a  class  of  naturally  occurring  polycrystal- 
line,  polyphase  ceramics.  Any  explanation  of  the  strength 
and  toughness  of  one  class  of  material  should  be  relevant 
to  the  other  class.  A  comparison  of  properties  taken  from 
the  literature  is  shown  in  Table  I.  Note  that  the  ceramics 
have  markedly  higher  unnotched  properties  while  the  notch 


References  are  given  on  page  55. 
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toughnesses,  R,  are  comparable.  In  fact,  if  average 
values  were  used  Instead  of  extreme  values  In  Table  I, 
the  rocks  would  be  found  to  be  tougher  than  the  ceramics, 
reversing  the  unnotched  behavior. 

The  large  values  of  fracture  energy  of  rock  are  of  particular  Interest. 

(8) 

Perkins  and  Bartlett  examined  this  problem  and  suggested  that  plastic  defor- 
matlon  occurs  at  the  tip  of  a  crack  In  rocks,  In  a  fashion  similar  to  metals. 

This  process  can  dissipate  large  amounts  of  energy  and,  If  It  occurs  In  rocks, 
can  account  for  the  difference  In  toughness  between  a  rock  and  Its  constituents. 

An  alternative  Interpretation  has  been  offered  by  Moavenzadeh  and  Kuguel^-  based 
on  their  work  on  concrete.  They  suggest  that  the  energy  dissipated  per  unit  crack 
area  formed  is  no  greater  In  rocks  than  in  Its  individual  constituents  but  only 
appears  greater  because  large  amounts  of  area  are  generated  by  microcracks  In 
regions  adjacent  to  the  main  fracture.  This  view  has  al^o  been  recently 
supported  by  Friedman,  et  al^11^  from  work  on  two  limestones,  sandstones,  and  a 
quartzite.  This  latter  interpretation  Is  also  strongly  supported  by  the  results 
presented  In  this  report. 

This  report  describes  research  aimed  at  overcoming  these  problems.  To 
this  end,  the  past  year's  work  aimed  at  a  detailed  description  of  the  failure  be¬ 
havior  of  two  sedimentary  rock  types:  Salem  limestone  and  Berea  sandstone.  The 
approach  conblned  the  measurement  of  R,  the  energy  dissipated  In  fracture,  with 
high  resolution  fractography  and  petrography  of  broken  rock  surfaces  and  the 
material  adjacent  to  the  fracture,  and  with  acoustic  detection  of  microcrack  form¬ 
ation. 


*  Also  called  micro fractures  and  en-echelon  microcracking. 
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EXPERIMENTAL  PROCEDURE 

The  energy  dissipation  rates  for  crack  extension  In  the  limestone 
and  sandstone  were  measured  by  means  of  the  double-cantilever-beam  (DCB)  test 
specimen  and  loading  arrangement  shown  In  Figure  1.  The  specimen  dimensions 
were  1.6  x  7.6  x  17.8  cm  with  a  slot  tipped  with  a  0.4-mm-wlde  sawcut  placed  to 
facilitate  starting  a  crack  having  the  desired  orientation  relative  to  the 
bedding  plane.  Prior  to  testing,  the  specimens  were  baked  at  110°C  In  order 
to  minimize  the  possible  effects  of  differences  In  moisture  content  from  one 
specimen  to  another.  A  few  specimens  were  tested  after  a  12-hour  immersion  In 
tap  water.  To  measure  the  growth  of  the  crack  a  parallel  array  of  conductive, 
lines  was  placed  on  each  specimen.  The  electrical  resistance  of  the  array 
was  monitored  during  testing  and,  as  the  crack  severed  the  lines,  the  corre¬ 
sponding  resistance  change  revealed  the  crack-tip  location.  Because  of  the 
extreme  brittleness  of  the  rock,  It  was  difficult  to  find  a  conductive  material 
sufficiently  brittle  so  that  the  strip  would  fell  reproduclbly  the  Instant  the 
crack  In  the  reck  passed  underneath.  Lines  drawn  on  the  surface  with  a  lead 
pencil  were  found  to  be  more  reliable  than  either  silver  paint  or  a  number  of 
vapor  deposited  metals,  although  some  uncertainty  In  the  crack  length  measure¬ 
ments  remained. 

Loads  were  applied  to  the  test  specimens  by  forcing  a  wedge  between 

a  pair  of  steel  pins  located  on  either  side  of  the  specimen.  Typical  test 

records  of  the  load  applied  to  the  wedge  and  the  resulting  crack  opening  uis- 

placerient  and  crack  growth  are  shown  In  Figure  2.  The  testing  arrangement 

favors  stable  crack  extension  when  the  starting  notch  is  a  sharp  crack,  and 

(121 

unstable  crack  propagation  when  the  starting  notch  is  blunted. 
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FIGURE  2.  A  TYPICAL  LOAD-DISPLACEMENT -RECORD  CRACK  LENGTH  FROM  A  WEDGE-LOADED  DCb  TEST.  In  tne 
initial  part  of  the  record  the  cracfc  has  not  yet  begun  to  grow  as  indicated  by  the 
linear  relation  between  the  load  and  displacement.  In  the  latter  part,  the  load  is 
gradually  decreasing  as  the  crack  extends  at  a  more  or  less  constant  rate. 
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The  fracture  energy  dissipation  rate,  R,  defined 
irreversibly  consumed  per  unit  area  of  crack  extension,  is 


as  the  energy 
derived  from  this 


test  by  way  of  the  following  energy  balance: 


R  '  dA  (1) 

where  dW  represent,  the  work  done  by  the  testing  machine  during  an  Increase 
of  the  crack  area  by  dA  and  dU  is  the  corresponding  change  In  strain  energy  In 
the  specimen.  Equation  (1)  applies  to  the  case  where  the  crack  Is  growing 
slowly  and,  therefore,  the  kinetic  energy  In  the  system  Is  negligible.  An 
analysis  based  on  Equation  (1)  has  been  developed  at  Battell.  which  allows  R 
to  be  computed  from  either  the  load-displacement  record  or  from  measured  dis¬ 
placements  and  crack  lengths.  Details  of  this  analysis  are  given  t„ 

Appendix  A.  By  convention,  R  Is  based  on  the  projected  surface  area: 

dA  -  bda 

Where  b  Is  the  specimen  thickness  and  a  la  the  crack  length.  As  shown  In  the 
next  section,  the  actual  crack  surface  Is  highly  Irregular,  and  the  true  surface 
area  Is  larger  than  assumed  In  calculating  R.  This  Is  an  Important  consider- 
ation  in  interpreting  the  results  of  these  tests. 

The  calculation  of  R  from  the  load-displacement  record  depends  on 
E,  Young  modulus.  The  value  of  E  employed  for  the  limestone  was  determined 
from  the  slope  (compliance)  of  the  load-displacement  record  for  a  DCB  specimen 
containing  a  crack  of  known  length  and  0,  the  modulus  independent  com¬ 
pliance  derived  from  elasticity  theory  (see  Equation  (A6)  in  Appendix 
A).  The  ratio  of  these  two  quantities  gives  the  value  of  E  as 
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1.6  x  10  dynes/cm  (2.3  x  10  psi) .  This  value  compares  quite  well  with 

the  zero  stress  value  obtained  from  uniaxial  compression  data  by  Hardy  and 
(13) 

Kim  on  similar  material.  The  elastic  modulus  applied  to  the  sandstone 

data  J.s  7  x  10*0  dynes/cm^  (1  x  10^  psi),  slightly  lower  than  the  geo-pressure 
(14) 

value  of  La  Mori  .  Because  the  sandstone  exhibits  a  greater  degree  of  non¬ 
linearity  (in  tension)  this  modulus  is  difficult  to  define  exactly.  These 
effects  and  their  implication  in  terms  of  a  fracture  model  for  these  materials 
are  discussed  later. 

Acoustic  emissions  accompanying  the  wedge  loading  were  monitored  in 
a  number  of  experiments.  The  emissions  were  detected  by  attaching  piezo¬ 
electric  crystals  to  the  specimen  surface.  The  crystals  were  acoustically 
coupled  to  the  specimen  by  means  of  a  viscous  resin.  The  crystals  have  a 
resonant  frequency  of  150  kHz  and,  therefore,  it  is  this  frequency  of  acoustic 
pulses  which  is  most  sensitively  detected.  In  tests  monitoring  acoustic 
emissions,  the  specimens  were  placed  on  a  rubber  pad  on  the  loading  platen 
and  teflon  sheets  were  used  between  the  wedge  and  pins.  These  precautions 
were  aimed  at  minimizing  the  stray  noise  from  specimen-platen  contact  or  from 
sliding  of  the  wedge  along  the  pin.  Tests  on  uncracked  specimens  showed  that 
these  precautions  were  effective,  and,  therefore,  acoustic  pulses  registered 
by  the  recording  equipment  during  an  actual  test  originate  in  the  bulk  of  the 
specimen,  presumably  in  the  region  of  high  stress  surrounding  the  slot  tip. 

The  crack  path  and  microstructural  features  of  the  fracture  surface 
were  examined  in  detail  in  both  materials.  These  examinations  involved  con¬ 
ventional  optical  (reflection)  microscopy,  scanning  electron  microscopy, 
transmission  electron  microscopy  of  Pt-C  shadowed  replicas,  and  optical 
(transmission)  microscopy  of  thin  slices. 


EXPERIMENTAL  RESULTS 


Load -Displacement  Crnck-Length  Measurements 

The  load  displacement  reco*  ds  for  both  the  limestone  and  the 
sandstone  can  be  divided  into  3  regions,  as  shown  in  Figure  2.  In  Region  I, 
the  load-displacement  record  is  more  or  less  linear  consistent  with  elastic 
behavior .  In  Region  II  the  record  then  becomes  increasingly  nonllnea.  as  a 
result  of  inelastic  behavior  in  the  region  of  high  stresses  surrounding  the 
notch. 

Typically,  no  gross  crack  extension  was  detected  in  Region  II, 
either  by  the  conductive  strips,  or  by  observation  with  a  traveling  micro¬ 
scope  at  20X  magnification.  Interior  sections  of  samples  loaded  into  Region 
II  also  failed  to  provide  evidence  of  gross  crack  extension.  Region  III  is 
characterized  by  decreases  In  load  with  increas.it  in  displacement,  as  is  expected 
for  a  growing  crack.  As  shown  in  Figure  2,  measurements  of  crack  extension  de- 
rived  from  the  conducting  strips  ,  show  that  the  onset  of  crack  extension  and 
Region  III  of  the  load  displacement  curve  approximately  coincide.  The  onset  of 
both  Region  II  and  Region  III  occurred  at.  lower  loads  on  the  sandstone  than  in 
the  limestone.  Furthermore,  and  contrary  to  expectations,  crack  growth  from 
the  relatively  blunt  starting  notches  (sawcuts)  proceeded  in  a  stable  rather 
than  In  an  unstable  fashion. 


*  These  were  confirmed  by  the  observations  with  the  traveling  microscope. 
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Fracture  Energy  Dissipation  Rates 

Values  of  the  fracture  energy,  R,  for  the  limestone  and  the  sandstone 
displayed  the  same  type  of  dependence  on  crack  length,  as  shown  in  Figure  3. 
From  the  first  evidence  of  crack  extension  the  R  values  increase  up  to  a  point 
and  then  level  out  at  a  relatively  constant  value  as  the  crack  continues  to 
grow.  This  level,  designated  as  the  asymptotic  limit  in  Table  II,  is  the  most 
■  mracteristic  value  of  R  for  the  rock  as  it  defines  the  R  associated  with  ex¬ 
tension  of  large  (relative  to  the  grain  size)  preexisting  flaws  in  the  rock.  A 
compilation  of  the  average  R  values  taken  from  the  upper  plateau  are  given  in 
Table  II.  For  the  limestone,  fracture  on  a  plane  parallel  to  the  bedding  plane 
(AY)  is  accompanied  by  the  lowest  energy  absorption  rate  while  cracking  on  a 
plane  perpendicular  to  the  bedding  plane  but  in  a  direction  parallel  to  it,  (BX) 
requires  substantially  more  energy.  The  energy  dissipation  rate  for  the  sand¬ 
stone  is  larger  than  for  the  limestone  and,  in  addition,  displays  a  different 
orientation  dependence.  The  easy  propagation  plsne  is  again  the  bedding  plane 
but  the  most  difficult  orientation  is  CZ,  the  plane  and  direction  perpendicular 
to  the  bedding  plane.  It  is  interesting,  but  probably  coincidental,  that  the 
ratios  of  R  for  the  easiest  and  most  difficult  orientations  are  roughly  the 

same  for  both  rocks.  R  values  for  comparable  rocks  obtained  from  other 
(8,11,15) 

sources  are  also  given  in  this  table.  While  they  did  not  report  the 

orientation  of  their  crack  plane  their  R  valors  are  comparable  to  toe  range  of 
data  from  this  study. 

A  few  tests  on  each  rock  were  conducted  in  1* juid  nitrogen  (-196BC). 

(9) 

Gilman  has  shown  that  at  this  temperature  calcite  undergoes  cleavage  with 
2 

an  R  ■  0.46  J/m  ,  close  to  the  specific  surface  energy;  consequently,  the 
energy  dissipated  by  plastic  deformation  is  negligible.  Since  the  crack  path 


Energy  Dissipation  Rate  , joules/m 
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FIGURE  3.  ENERGY  DISSIPATION  RATES  MEASURED  OVER  AN  INTERVAL  OF  CRACK  LENGTH 
FOR  THREE  ORIENTATIONS  OF  THE  CRACK  PLANE  AND  PROPAGATION  DIRECTION 
IN  SALEM  LIMESTONE.  The  letters  A,  B,  and  C  denote  orientation  of 
the  crack  plane  with  respect  to  the  bedding  plane  (plane  A)  and  X, 
Y,  and  Z  denote  crack  propagation  direction.  a/aQ  is  the  ratio  of 
crack  length  to  the  starting  slot  length. 
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TfBLE  II.  COMPILATION  OP  ENERGY  DISSIPATION  RATE  RESULTS 


Material 

Orientation^ 

Asymptotic- 
R  (joules/tn  ) 

Salem  Limestone 

This  Study 

dry(b)  (22°C) 

AY 

50-75 

BX 

195-230 

CZ 

90-110 

(c)  (-196»C> 

CZ 

85 

wet' 

CZ 

~  40 

Perkins  and  Bartlett ^ 
Friedman,  et  al(e' 

84 

55 

Berea  Sandstone 

This  Study 
dry<h)  (22°C) 

AY 

465-520 

BX 

600-750 

CZ 

1120-1580 

(-196°C) 

wet'0' 

AY 

BX 

150-270 
~  300 

Perkins  and  Krech^ 
Friedman,  et  al(e) 

176 

139 

(a)  Refer  to  Figure  3. 

(b)  After  baking  at  110°C  for  1  hour,  specimens  were  stored  in  a  dessicacor 

prior  to  testing.  * 

(c)  Specimens  were  soaked  in  tap  water  for  approximately  10  hours  prior  to 

L.wO  Llug  * 

(d)  Reference  8.  A  side  grooved  DCB  specimen  was  used  having  a  width  of 
15.2  cm  and  net  thickness  of  3.8  cm.  This  R  value  is  twice  the  value 
of  surface  energy  quoted  in  References  1,  3,  and  6  since  2V  «  R 
Orientation  not  reported. 

(e)  Reference  11 .  A  notch  bend  specimen  was  used  with  a  thickness  of 
1.4  cm  and  beam  height  of  2.5  cm.  Their  sandstone  is  listed  as 
Tennessee  sandstone.  These  values  were  derived  from  the  data  of 
ir-iedman,  et  al,  using  the  projected  surface  area.  Orientation 
not  reported. 

Reference  15.  Same  specimen  design  and  dimensions  as  in  Reference  1. 
'heir  sandstone  is  listed  as  Tennessee  sandstone.  Orientation  not 
reported. 
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in  both  the  limestone  and  sandstone  is  through  calcite,  a  comparison  of  the 
R  value  at  -196°C  with  the  room  temperature  value  should  indicate  the  degree 
to  which  plastic  deformation  contributes  to  the  room  temperature  fracture 
energy  dissipation  processes.  Th  results,  given  in  Table  II,  show  that  the 
R -values  for  limestone  and  sandstoe  are  still  ~  100X  to  1000X  the  value  obtained 
by  Gilman.  Thus,  plastic  deformation  cannot  account  for  the  magnitude  of  the 
R -values  obtained  either  at  room  temperature  or  at  -196°C.  The  authors  favor 
the  alternative  explanation  that  the  temperature  dependence  of  the  R  in  the 
sandstone  has  its  origin  in  the  difference  in  thermal  expansion  coefficients 
between  quartz  and  calcite.  On  cooling  the  sandstone  to  -196,'C,  the  difference 
in  thermal  expansion  between  the  two  components  in  this  rock  should  produce 
internal  stresses  which  may,  in  turn,  aid  the  applied  stresses  in  producing 
fracture.  The  limestone,  being  monomlnerallc ,  would  be  expected  to  undergo 
much  less,  if  any,  differential  thermal  expansions  leading  to  internal  stresses. 

The  introduction  of  water  appears  to  produce  a  moderate  decrease 
in  R  consistent  for  both  rocks.  The  fact  that  the  R  values  for  both  rocks 
decrease  by  the  same  factor  is  probably  related  to  the  fact  that  the  crack 
path  is  through  calcite  in  each  case.  Hence,  environmental  effects  could  be 
expected  to  produce  similar  changes  in  the  fracture  behavior  of  these  two 
materials.  There  is,  however,  ecme  uncertainty  in  the  magnitude  of  these  R 
values  for  the  wet  condition  because  of  difficulty  encountered  with  the 
conductive  strips  used  as  crack  length  indicators.  Until  a  more  reliable 
procedure  can  be  found  for  detecting  the  position  of  the  crack  in  the  wet 
specimens  these  values  should  serve  only  as  reasonable  estimates. 
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Acoustic  Emission 

The  tests  in  which  acoustic  emissions  were  monitored  revealed  some 
interesting  information  about  the  fracture  mode  In  these  materials.  Examples 
of  load-time  curves  superposed  on  an  acoustic  activity  versus  time  bar  graph 
are  given  in  Figures  4  and  5  for  the  limestone  and  sandstone,  respectively. 
Significant  acoustic  activity,  which  quite  likely  derives  from  cracking,  is 
detected  at  loads  well  below  the  maximum  load  which  corresponds  to  the  point 
of  optically  observable  crack  extension.  Furthermore,  the  first  activity 
above  background  occurs  at  about  the  load  levels  at  which  the  load-displacement 
record  becomes  nonlinear.  This  evidence  coupled  with  other  results  described 
below  points  to  the  development  of  numerous  cracks  in  the  volume  surrounding 
the  slot  tip  which  may,  as  a  result,  account  for  the  gradual  decrease  in  slope 
in  the  early  stages  of  the  test.  In  figure  4  the  Advance  of  the  crosshead  was 
stopped  several  times  during  this  test.  Upon  stopping  the  crosshead,  the 
acoustic  activity  did  not  terminate  immeddately,  but  instead  the  count  rate 
decayed  gradually,  as  did  the  load,  indicating  some  additional  growth  of  the 
crack.  This  behavior  of  crack  growth  under  fixed  displacement  conditions  is 
often  observed  in  metals  and  is  referred  to  as  subcritical  cracking.  It  is 
indicative  of  a  reaction  between  the  environment  and  the  exposed  surface  at 
or  very  near  the  tip  of  the  crack. 

During  the  periods  of  crosshead  advance  the  acoustic  activity  of  the 
limestone  was  similar  in  level  to  that  of  the  sandstone.  While  thitf  may  imply 
that  the  amount  of  cracking  is  the  same  in  both  rocks,  it  would  be  incorrect 

to  conclude  that  this  is  necessarily  true  because  differences  in  acoustic  wave 
generation  and  transmission  properties  between  the  two  rocks  can  markedly  affect 
the  amplitude  of  the  signals  received  and,  hence,  the  count  rate. 


Load  .Sbs 
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FIGURE  4.  AN  EXAMPLE  OF  THE  RATE  OF  ACOUSTIC  PULSES  RECEIVED 
DURING  TESTING  OF  A  SALEM  LIMESTONE  SPECIMEN.  The 
count  rate  is  correlated  with  the  wedging  load 
applied  to  the  specimen.  The  noise  or  background 
level  is  approximately  30  counts/minute.  The  cross¬ 
head  advance  of  the  testing  machine  was  stopped  at 
the  time  designated  by  the  double-ended  arrows  and 
restarted  at  the  single-ended  arrows. 


Count  Rate,  countsArunule 
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FIGURE  5.  AN  EXAMPLE  OF  THE  RATE  OF  ACOUSTIC  PULSES  RECEIVED  DURING 

TESTING  OF  A  BEREA  SANDSTONE  SPECIMEN.  Here  the  background 
level  is  about  25  count s /minute  and  the  crosshead  was  not 
stopped  during  the  test. 
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Microscopic  Studies  of  the  Salem  Limestone 

The  microptructure  of  the  Salem  limestone  specimens  consists  of  an 
aggregate  of  particles  Identifiable  as  calclferous  fossils  or  fragments  of 
fossils  together  with  'oolites.  The  particles  are  loosely  packed  and  consider¬ 
able  pore  space  Is  observable.  The  diameters  of  the  particles  in  this  aggregate 
are  principally  in  the  range  of  0.5  mm  to  2.0  mm.  At  low  magnification  (~  30X) 
the  particles  appear  to  be  cemented  by  a  white  powdery  material  (the  matrix). 

At  higher  magnifications  the  granular  character  of  the  matrix  is  observable  as 
can  be  seen  in  Figure  6  which  shows  electron  micrographs  of  a  replica  made  from 
a  polished  and  lightly  etched  surface.  Both  views  in  Figure  6  show  the  structure 
within  the  matrix  or  at  the  interface  of  a  larger  particle  (oolite  or  fossil 
fragment).  The  fine  particles  comprising  the  matrix  appear  to  have  a  diameter 
of  about  1  fi.  The  reflecting  optical  micrographs  in  Figure  7  display  the 
heterogeneous  nature  of  the  microstructure  and  also  show  the  path  of  the  crack. 
The  micrographs  in  Figure  7  were  obtained  from  a  polished  surface  of  a  specimen 
with  a  crack  of  the  AY  orientation.  The  specimen  was  first  prepared  for  obser¬ 
vation  by  infiltrating  the  crack  and  the  interconnecting  pores  with  an  epoxy 
in  order  to  prevent  tearing  out  of  loose  material  during  polishing.  These 
micrographs  indicate  that  the  crack  propagates  most  often  through  the  matrix 
described  above.  There  is,  however,  considerable  fracture  of  the  fossil  debris 
and  oolites  as  can  be  seen  in  Figure  7.  Note  also  that  because  the  crack  pre¬ 
dominantly  follows  the  matrix,  the  path  is  highly  irregular  and,  therefore,  the 
actual  surface  area  created  during  fracturing  is  larger  than  the  area  referred 
to  in  Equation  (2)  and  on  which  the  calculation  of  R  is  based. 


V  fjr 
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Proceeding  to  higher  magnifications  of  the  fracture  surface  repeals 
more  clearly  the  highly  irregular  character  of  the  surface  topography.  Figures 
8a,  b,  and  c  are  a  sequence  of  stereo  pairs*  of  progressively  increasing  magni¬ 
fication  made  by  scanning  electron  microscopy  of  an  area  of  the  fracture  surface 
of  a  Salem  limestone  specimen.  In  Figure  8a,  the  outline  of  a  number  of  milli¬ 
meter  size  particles  are  observable.  Figure  8b  reveals  an  area  where  the  crack 
has  passed  partially  through  the  matrix  material  as  well  as  through  larger 
particles.  The  fracture  surfaces  of  the  larger  particles  (whether  fossil 
fragments,  precipitated  crystals,  or  oolites)  display  distinct  patterns, 
reminiscent  of  cleavage  of  single  crystals.  Verification  that  these  are  single 
crystals  of  calcite  or  at  least  highly  oriented  polycrystals  is  discussed  in 
reference  to  Figure  9.  Ac  high  magnification  the  fine  micron  size  particles 
comprising  the  matrix  appear  to  be  faceted  as  shown  in  Figure  8c  which  suggests 
that  these  also  may  be  single  crystals.  The  utility  of  the  scanning  electron 
microscope  derives  f i  om  the  large  depth  of  focus  it  prrvides  at  all  levels  of 
magnification.  A  good  example  of  this  capability  is  shown  in  Figure  8d  which 
reveals  a  pore  broken  open  during  fracture  and  which  contains  a  lining  of 
nearly  perfect  crystals  of  calcite  or  perhaps  magnesite. 

In  order  to  resolve  features  as  small  as  20  angstroms,  plastic 
replicas  were  made  of  the  fracture  surfaces.  A  platinum-carbon,  mixture  was 
then  vapor  deposited  onto  the  plastic  replica  which,  in  turn,  was  then  dis¬ 
solved  in  acetone ,  The  Pt-C  films  are  placed  in  an  electron  microscope  in 
which  an  electron  beam  passes  through  the  films  and  the  resuiting  image  can  be 
displayed  on  a  phosphor  screen  or  pho-graphed .  During  the  replication  process 

*  A  three-dimensional  effect  is  obtained  when  viewing  these  micrograph  pairs 
through  a  stereo  viewer.  If  such  a  veiwer  is  available  the  pages  should  be 
removed  and  slit  down  the  center  in  order  to  adjust  the  spacing  between  the 
micrographs  for  the  viewer  used. 
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STEREO  PAIRS  OBTAINED  BY  SCANNING  ELECTRON  MICROSCOPY  OF  FRACTURE 
SURFACE  IN  SALEM  LIMESTONE  (Continued). 
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FIGURE  8.  STEREO  FAIRS  OBTAINED  BY  SCANNING  ELECTRON  MICROSCOPY  OF  FRACTURE 
SURFACE  IN  SALEM  LIMESTONE  (Continued) . 
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of  the  limestone,  small  fragments  on  the  fracture  surface  were  removed  and 
could  be  observed  in  the  final  replica.  For  example,  Figures  9a  and  9b  show  the 
outlines  of  numeorus  micron  size  particles  comprising  the  matrix,  as  well  as 
the  impression  in  the  replica  of  areas  of  the  surface.  Contrast  the  irregular 
character  of  this  surface  with  the  relatively  flat  cleavage  surface  of  one  of 
the  millimeter  size  particles  in  Figure  9c.  The  distinct  feature  in.  Figure  Sc 
is  known  as  a  river  pattern  and  occurs  when  the  crack  passes  from  one  single 
crystal  into  another  with  a  very  slight  difference  in  crystallographic  orien¬ 
tation.  It  was  also  discovered  that  thin  flakes  of  material  were  removed  from 
these  cleavage  surfaces.  These  flakes  were  thin  enough  to  be  transparent  to  the 
electron  beam  thereby  relevating  the  internal  structure  as  can  be  seen,  for 
example,  at  ultra-high  magnification  (108,000X)  in  Figure  9d.  In  addition, 
electron  diffraction  patterns  can  be  obtained  from  these  flakes,  an  example  of 
which  is  shown  in  Figure  9e.  These  patterns  can  be  used  to  make  relatively 
precise  determinations  of  crystal  structure.  The  pattern  in  Figure  9e  has  been 
analyzed  and  was  found  to  correspond  with  thf  structure  of  single  crystal  calcite. 


Microscopic  Studies  of  the  Berea  Sandstone 


j.he  microstructure  consists  of  remarkably  uniform-sized  quartz  grains 
(0.3  mm  dia.),  cemented  by  calcite.  A  considerable  number  of  voids  are  evident 
and  Hand  in ,  et  al(16^  report  a  volume  fraction  of  pores  of  187..  The  main  crack 
propagates  through  the  interconnecting  calcite  layers.  Figure  10  provides  some 
evidence  that  prior  to  and  accompanying  the  propagation  of  the  main  crack  there 
is  extensive  cracking  occurring  in  a  relatively  large  volume  of  material 


adjacent  tc  the  crack.  During  vacuum  impregnation  of  the  sandstone,  with  epoxy 
resin,  penetration  was  typically  only  2  to  3  mm  into  the  surface.  However, 


FIGURE  9.  HIGH  RESOLUTION  ELECTRON  MICROSCOPE  FRACTOGRAPHS  OBTAINED  FROM  Pt-C  SHADOWED 
REPLICAS  OF  THE  FRACTURE  SURFACE  OF  SALEM  LIMESTONE. 
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FIGURE  9.  HIGH  RESOLUTION  ELECTRON  MICROSCOPE  FRACTOGRAPHS  OBTAINED  FROM  Pt-C  SHADOWED 
REPLICAS  OF  THE  FRACTURE  SURFACE  OF  SALEM  LIMESTONE  (Continued). 


more  permeable  as  can 


around  the  tip  shown  in  Figure  9  the  material  is  much 

be  seen  by  the  boundary  of  the  .poxy  extending  farther  than  1  cm  from  the  slot 

tip.  Furthermore,  during  sectioning  of  specimens  containing  a  crack,  it  was 

noticed  that  quartz  grains  adjacent  to  the  crack  were  quite  loose  and  could  be 
easily  removed. 

On  the  fracture  surface  at  relatively  low  magnifications,  individual 
quartz  grains  are  readily  observable  as  can  be  seen  on  the  stereo  pairs  ob¬ 
tained  by  scanning  electron  microscopy  in  Figure  11.  The  grains,  exhibit  both 
irregular  surfaces  and  planar  facets.  There  is  also  some  indication  here  of 
very  narrow  cracks  in  the  seams  between  individual  grains.  In  fact,  at  100X 
(Figure  11a)  the  SEM  fractographs  give  an  &r,pes  ;..<nce  of  loosely  adherent 
rubble.  This  appearance  is  somewhat  consistent  with  the  fact  that  after  the 
sandstone  specimens  were  completely  broken  open  a  large  number  of  particles 
would  fall  from  the  fracture  si cfaces  and  many  more  could  be  dislodged  with  a  slight 
tap.  Figure  11  b  shows  the  only  evidence  of  cracking  of  the  quartz  grains  that 
was  found.  This  shows  an  area  of  conchoidal  fracture  with  a  small  chip  partially 
broken  away  from  the  underlying  quartz  grain.  In  Figure  11c  the  calcite  cement 
can  be  identified  as  the  aggregate  of  small  rather  blocky  particles.  Another 
area  of  calcite  is  shown  at  higher  magnification  (2000X)  in  Figure  lid.  The 
calcite  appears  as  a  well-defined  polycrystalline  aggregate  displaying  many 
surfaces  identifiable  as  cleavage  surfaces.  These  cleavage  surfaces  are  dis¬ 
tinguished  by  a  river  pattern  on  the  surface,  an  excellent  example  of  which  is 
shown  in  Figure  12a  which  is  obtained  by  electron  microscopy  of  a  replir*  f  om 
the  fracture  surface.  It  cannot  be  firmly  established  whether  cleavage  areas 
in  Figure  12a  are  in  the  quartz  or  the  calcite  but  their  general  appearance  is 
similar  to  the  cleavage  morphology  found  in  tie  limestone,  and  it  is  likely  that 
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FIGURE  II.  STEREO  PAIRS  OBTAINED  BY  SCANNING  ELECTRON  MICROSCOPY  OF  THE  FRACTURE  SURFACE  OF  BEREA  SANDSTONE. 


FIGURE  11.  (Continued) 


■  -  1  V  j 

\  >  +m£L*zI^3£7i 

'  .  y^;7  i‘^n  *  1 

j'  AiL  v  !S 

L 

FIGURE  12.  HIGH  RESOLUTION  ELECTRON  MICROSCOPE  FRACTOGRAPHS  OBTAINED  FROM  Pt-C  SHADOWED  REPLICAS 
OF  THE  FRACTURE  SURFACE  OF  BEREA  SANDSTONE:  (a)  areas  of  distinct  cleavage  and  (b) 
secondary  silica  deposits. 
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these  cleavage  areas,  which  were  cc  mon  on  the  replicas,  predominantly  represent 
the  fracture  surface  of  the  calcite  cement.  The  other  feature  commonlv  observed 
is  shown  in  Figure  12b  and  probably  is  representative  of  the  surface  of  quartz 
grains.  These  blocky  crystals  exhibit  typical  habit  planes  of  quartz  and  are 
similar  to  features  found  on  the  surface  of  quartz  grains  by  Kirkpatrick  and 

(17) 

Summerson  .  Therefore,  these  surfaces  probably  form  the  pore  walls  and  are 
not  produced  as  the  result  of  fracture. 

When  examining  the  stereo  pairs  one  obtains  a  very  vivid  portrayal 
of  highly  uneven  relief  of  the  fracture  surface  in  these  materials.  Some  quanti¬ 
tative  measurements  of  this  relief  were  made  from  Figures  11a,  c,  and  d  by  means 
of  a  Hilger  and  Watts  stereoscope.  The  resulting  elevation  profiles  are  shown 
in  Figures  13a,  b,  and  c  along  with  an  inset  which  indicates  che  position  of  the 
profile  on  each  of  the  SEM  fractographs .  At  100X,  Figure  13a,  the  profiles  of 
individual  quartz  grains  are  easily  recognized  and  going  to  500X  offers  little 
additional  information  concerning  the  topography.  However,  at  2000X  in  Figure 
13c,  the  faceted  character  of  the  calcite  cement  becomes  obvious,  a  feature 
which  cannot  be  distinguished  in  the  profiles  at  lower  magnification. 

The  great  usefulness  of  these  elevation  profiles  is  that  they  provide 
a  means  for  estimating  the  exposed  curface  area  by  simply  squaring  the  actual 
length  of  the  line  representing  the  profile.  The  overall  profile  length  can 
be  most  adequately  obtained  from  the  magnifications  at  10X,  100X,  and  2000X. 

The  profile  length  is  longer  than  the  horizontal  distance  by  a  factor  of  approxi¬ 
mately  1.2,  1.6,  and  1.7  at  each  of  these  magnifications,  respectively.  There¬ 
fore,  the  total  profile  length  determinable  up  to  a  magnification  of  2000X  is 
simply  the  product  of  these  factors  or  about  3.3  times  the  horizontal  distance. 
This  gives  an  exposed  surface  area  which  is  about  11  times  the  projected  surface 


area . 


Relative  Height, 
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(a) 

FIGURE  13.  ELEVATION  PROFILES  OF  TTE  FRACTURE  SURFACE  OF  BEREA  SANDSTONE. 

These  profiles  were  measured  using  a  Hilger  and  Watts  stereo¬ 
scope  from  the  stereo  pairs  in  Figure  11.  The  inset  shows 
the  location  of  each  profile  on  the  corresponding  figure. 

(a)  From  Figure  11a  at  100X. 

(b)  From  Figure  11c  at  500X. 

(c)  From  Figure  lid  at  2000X. 
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DISCUSSION 

The  results  obtained  in  this  study  favor  the  Interpretation  of 
Moavenzadeh  and  Kuguel^10^  and  Friedman  and  coworkers^11^,  which  ascribes  the 
relatively  large  fracture  energies  of  rocks  to  microcracks  generated  in  large 
numbers  throughout  tne  highly  stressed  region  surrounding  the  crack  tip.  The 
two  alternative  energy  sinks:  the  main  fracture  surface  and  plastic  formation 
derive  little  support  from  this  study.  The  detailed  measurements  of  the  sur¬ 
face  topography  at  magnifications  up  to  2000X  indicate  that  the  total  surface 

area  is  about  llx  the  projected  area  of  fracture.  Even  though  this  factor  is 

* 

larger  than  values  of  1.5  tc  3  and  1.8  to  4  that  have  been  measured  by 

Friedman,  et  al^11^,  for  Severn!  rocks  and  estimated  by  Wiederhon/18^  from 

geometrical  considerations,  it  is  not  large  enough  to  account  lor  the  fracture 

energy.  Note  that  the  measured  cleavage  fracture  energy  of  calcite  at  -196°C 
.  „  . .  .  2(9) 

is  0.46  j/m  ;  11  times  this  value  is  still  only  ~  37.  of  the  R-value  of  the 
sandstone  at  this  temperature.  There  is  also  little  basis  for  believing  that 
the  small  crystallites  in  the  calcite  cement  of  the  limestone  and  the  sandstone 
are  more  likely  to  dissipate  fracture  energy  by  undergoing  plastic  deformation 
at  -196°C  than  the  massive  crystal  tested  by  Gilman.  For  one  thing,  the  iso¬ 
lated  crystals  extracted  from  the  fracture  produced  at  room  temperature  dis¬ 
played  no  dislocations  when  viewed  in  transmission.  The  fact  that  the  appear¬ 
ance  of  the  fracture  calcite  cement  is  exceedingly  granular  in  the  limestone 
(see  Figure  8)  implies  that  the  crystallite  boundaries  are  weaker  than  the 
cleavage  plane,  that  dislocation  generation  is,  if  anything,  less  likely,  and 
that  the  fracture  energy  (per  unit  surface  area)  could  be  somewhat  smaller  than 
the  value  measured  by  Gilman.  It  should  also  be  noted  that  the  contribution  of 

*  These  numbers  are  based  on  examinations  at  100X  magnification  and  may  be 
compared  with  the  value  3.7  obtained  in  this  study  at  100 -magnificat ion. 
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plastic  deformation  to  the  frac^re  energy  of  polycrystalline  MgO,  a  so-called 
ductile  ceramic  is  only  18  J/m2  ,  and  would  not  account  for  the  fracture 
energies  of  the  rock  samples.  It  thus  appears  that  plastic  deformation  is  not 
a  major  energy  sink,  either  at  room  temperature,  or  at  -196°C.  The  substan¬ 
tial  reduction  of  the  R-value  for  the  sandstone  with  temperature  can  be 
accounted  for  in  terms  of  internal  stresses  arising  from  differential  thermal 
expansion  and  microcracking  as  discussed  in  the  next  section. 

In  contrast,  microcrack  formation  is  supported  by  a  number  of  obser¬ 
vations  of  the  present  study,  as  well  as  the  work  of  Moavenzadeh  and  Kugel(10) 
and  Friedman  and  coworkers^^s 

^coystic  Emissions.  Large  amounts  of  acoustic  emissions  are 
detected  in  Region  II,  before  the  extension  of  t\e  main  crack 
is  observed.  It  is  difficult  to  account  for  these  emissions 
except  as  instances  of  microcracking  in  the  region  surrounding 
the  initial  slot  tip. 

^  Nonlinear  Compliance.  It  was  generally  found  that  acoust  c 
activity  initiated  at  about  the  beginning  of  nonlinearity  in 
the  load -displacement  record.  The  correlation  is  consistent 
since  the  formation  of  numerous  microcracks  will  reduce  the 
effective  elastic  modulus  of  the  specimen. 

(3)  Permeability.  The  increased  permeability  to  the  epoxy  resin  in 
the  sandstone  is  another  indicator  of  extensive  cracking  adjacent 
to  the  initial  slop  tip  and  main  fracture  path.  Increases  in  per¬ 
meation  rate  by  a  factor  of  two  to  three  are  suggested  by  the  resin 
penetration  distances  shown  in  Figure  10.  It  should  be  noted  that 
this  specimen  was  impregnated  under  zero  load  and  some  microcracks 
will  be  partially  or  completely  closed,  therefore,  this  permea¬ 
bility  Increase  is  lower  than  might  be  expected  under  load. 
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(4)  X-Rav  Studies.  Friedman  and  coworkers demonstrate  a  relaxation 
of  residual  elastic  strains  both  in  the  quartz  grains  of  Tennessee 
sandstone  in  a  zone  as  deep  as  5  mm  adjacent  to  the  fracture  surface 
of  broken  bend  specimens  and  in  the  quartzite.  Such  relaxations  are 
believed  to  be  the  result  of  rupture  of  the  interconnecting  calclte 
cement. 

(5)  Direct  Observation.  Quantitative  measurements  by  optical  micro¬ 
scopy  by  Moavenzadeh  and  Kugel  show  that  the  crack  area  contains 
numerous  microcracks  adjacent  to  the  main  crack  in  broken  concrete 
specimens  equivalent  to  as  much  area  ac  20  times  the  cross-sectional 
area  of  their  specimens.  Microcracking  at  the  grain  boundaries  of 
Chichowee  quartzite  has  also  been  revealed  by  Friedman  and  coworkers 
with  a  silver  nitrate  staining  technique. 


Additional  evidence  that  the  extent  and  frequency  of  Che  microcracking 
is  consistent  with  the  observed  R-values  is  derived  from  the  compliance  measure¬ 
ments.  This  it.  done  by  assuming  that  the  point  of  deviation  from  linearity  in 
the  load-displacement  record  is  coincident  with  the  onset  of  extensive  micro- 
cracking,  and  by  estimating  the  total  area  of  microcracks  on  this  basis.  From 

a  .alysis  developed  in  Appendix  B  we  estimate  the  total  surface  area  associated 

* 

with  freshly  produced  microcracks  in  a  DCB  specimen  as 


where 


A 


4rra 


.3(5-20).  +  , 

2S  18(S-2a)  1 


(3) 


*  Equation  (3)  does  not  contain  the  area  of  preexisting  voids  and  cracks. 
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a  "  slot  length  in  the  specimen 

h  ■  beam  height  (one-half  total  svidth) 

b  ■  specimen  thickness 

a  °  average  area  per  microcrack 

ac  "  average  critical  tensile  stress  required  to  crack  the 

calcite 

F  ■  applied  load. 

Equation  (3)  la  obtained  by  approbating  the  DCB  apeclmen  by  two 
cantilever  beans  anchored  to  a  rigid  foundation.  This  approximation  la  quite 
adequate  In  view  of  two  other  assumptions,  namely:  (.)  the  mlcrocrack.  are 
noninteracting  and  (b)  ffc  la  the  maximum  tensile  .tress  which  can  be  supported 
by  the  rock.  Denoting  the  bracketed  term  on  the  right  aide  of  Equation  (3) 

F,  average  values  of  ?  obtained  at  maximum  load  In  teats  of  three  llme.tone 
specimens  and  of  six  sandstone  specimens  are  given  In  Table  III. 


TABLE  III,  ESTIMATED  TOTAL  MICROCRACK  SURFACE  AREAS 


Rock 

Estimated  Single 
Microcrack2 
Area  (a)  -m 

Average  F 
at  Max. 
Load-in 

Total  Microcrack  Area  at 
Max.  Load  from  Eq .  3 

m2 

Limestone 

10'6 

3.0  x  10 “3 

0.154 

Sandstone 

10  “7 

5.4  x  10 "3 

2.76 

Ihe  "lcr°cr*ck  slze  m  estimated  on  the  basis  that,  when  first  / 
formed,  the  .res  of  each  microcrack  will  be  approximately  equal  to  the  particle 
cross-sectional  area.  The  results  Indicts  that  st  maximum  load,  when  gross 
crack  extension  Initiates,  the  accumulated  are.  of  freshly  forced  microcracks 


2  2 

in  the  limestone  is  about  0.15  m  and  about  2.8  m  in  the  sandstone.  While 
these  calculated  areas  should  serve  only  as  rough  estimates  they  certainly 
indicate  that  the  amount  of  microcrack  area,  as  computed  simply  from  the 
loss  of  compliance  of  the  specimen,  is  truly  very  large.  By  virtue  of  the 
inverse  dependence  of  A  on  a,  Equation  (3)  also  predicts  that  microcracking 
should  be  more  extensive  in  the  sandstone  than  in  the  limestone  by  about  a 
factor  of  10-20.  The  measured  R-values  for  the  sandstone  are  about  2  to  8 
times  larger  than  for  the  limestone  whicn  is  in  very  reasonable  agreement 
with  the  prediction  of  Equation  (3)  considering  the  fact  that  <7  may  be 
somewhat  different  in  these  two  rocks.  Thus,  while  the  model  on  which 
Equation  (3)  is  based  is  crude  (and  could  certainly  be  improved),  it  does 
give  a  reasonable  account  of  the  observed  behavior.  Furthermore,  it  esta¬ 
blishes  an  important  link  between  microstructure  end  the  fracture  behavior 
of  rocks. 

The  conceptual  picture  describing  the  development  of  microcracks 
at  the  tip  of  a  machined  slot  in  these  two  rock  materials  is  summarized 
schematically  in  Figure  14.  Point  (a)  corresponds  to  the  unloaded  conditions. 
A  few  fresh  microcracks  introduced  during  cutting  surround  the  slot.  As  the 
load  is  applied,  some  cracking  occurs  on  the  weakest  planes  very  near  the  slot 
tip  as  in  (b) .  At  this  point  too  few  microcracks  have  formed  to  produce  a 
detectable  nonlinearity  in  the  load-displacement  record  and  they  formed  at 
stresses  too  small  to  produce  an  acoustic  signal  of  detectable  amplitude. 

With  further  increases  in  load  the  microcracked  region  becomes  larger  and 
stronger  interfaces  and  linkages  rupture  producing  a  much  more  intensely 
cracked  volume  surrounding  the  slot  tip  as  envisaged  in  (c)  .  The  loading 
record  is  now  nonlinear  (Region  II)  and  the  rupture  of  the  interfaces  at  these 
load  slot  tip  stress  levels  produces  detectable  acoustic  signals.  There  may 
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be  locations  along  the  leading  edge  of  the  slot  where  some  microcracks  have 
linked  together  and  have  in  turn  joined  the  slot  tip.  However,  it  seems 
more  appropriate  to  simply  refer  to  this  region  as  a  zone  of  damage  which 
is  still  capable  of  supporting  load  and,  therefore,  the  conditions  for  ex¬ 
tension  of  the  main  crack  are  still  unsatisfied.  Further  increases  in  load 
extends  the  damage  zone  somewhat,  but  at  the  slot  tip  rupture  of  the  remain¬ 
ing  strongest  linkages  occurs  and  gross  crack  extension  begins  at  about  the 
point  of  maximum  load  (onset  of  Region  III) .  At  this  stage  expansion  of  the 
damage  zone  may  still  be  more  rapid  than  the  growth  of  the  main  ^rack.  This 
is  suggested  by  the  observation  that  R  increases  in  the  initial  stages  of 
crack  extension.  Finally,  however,  the  rate  of  propagation  of  the  main 
crack  tip  becomes  equal  to  the  rate  of  growth  of  the  damage  zone  extending 
ahead  of  the  crack  tip.  Thus,  a  steady-state  situation  is  achieved  in  which 
the  main  crack  extends,  pushing  the  damage  zone  ahead  of  it  at  the  same  rate 
as  is  depicted  in  (d) .  R  becomes  equal  to  the  asymptotic  limit  determined 
by  the  size  of  the  damage  region  and  the  rate  of  microcrack  production  within 
it.  Inherent  in  this  model  is  the  supposition  that  the  various  interfaces 
which  act  as  potential  microcrack  sites  exhibit  a  broad  range  of  load-carrying 
capability.  If  each  of  these  interfaces  could  support  the  same  tensile  load, 
microcracking  would  be  confined  to  the  slot  tip  and  the  development  of  a  large, 
energy -absorbing  damage  zone  prior  to  the  growth  of  the  main  crack  would  not 
be  possible.  It  is  interesting  to  note  that  similar  concepts  and  the  substi¬ 
tution  of  crazed  regions  for  microcracks  have  recently  been  elaborated  on  by 
(20  21) 

Hull,  et  al  ’  ,  to  explain  the  toughness  of  glassy  polymers. 

As  point"!  out  by  Perkins  and  Bartlett^  one  commonly  finds  a  large 
discrepancy  between  the  measured  tensile  strength  of  rocks  and  the  theoretical 
tensile  strength  based  on  a  simple  Griffith  theory  and  measured  R-values.  As 
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an  example,  if  the  rock  contains  cracks  of  length  (or  diameter  for  periny- 
diaped  cracks),  2a,  oriented  perpendicular  to  an  applied  uniaxial  stress  then 
the  tensile  strength  is 


where  P  is  a  geometry  dependent  constant  of  order  unity.  Table  IV  compares 
measured  tensile  strengths  with  calculated  values  based  on  R-measurements  on 
the  same  materials.  It  shows  that  measured  values  of  strength  fall  far  below 
the  predicted  levels . 

The  conceptual  picture  of  the  zone  of  damage  described  above  offers 
a  possible  explanation  for  this  discrepancy.  This  zone,  while  providing  a 
neans  of  dissipating  large  amounts  of  energy,  also  serves  another  very  impor¬ 
tant  function.  The  formation  of  numerous  microcracks  in  this  zone  diminishes 
the  stress  concentration  at  the  tip  of  the  main  crack  or  slot.  The  consequences 
of  this  process  are  precisely  the  same  as  in  a  metal  in  which  a  plastic  zone 
developes  in  response  to  large  stresses  at  the  crack  tip.  A  typical  relation 
between  fracture  stress  and  crack  size  for  a  metal  is  represented  in  Figure 
15.  For  long  crack  lengths  the  fracture  stress  varies  with  crack  length  in  a 
manner  described  by  Equation  (4)  since  the  R  associated  with  onset  of  rapid 
crack  extension  is  a  material  property  and,  hence,  constant.  However, 

Equation  (4)  only  applies  as  long  ns  the  size  of  the  plastic  zone  is  small 
relative  to  the  dimensions  of  the  structure.  As  the  crack  size  is  reduced 
the  fracture  stress  increases  until  it  approaches  the  yield  strength  and  wide¬ 
spread  plastic  deforr^tion  of  the  structure  occurs  prior  to  fracture. 

Equation  (4),  wh.'.ch  is  based  on  a  linear  elastic  analysis,  no  longer  applies 
and  fracture  stress  levels  out  and  becomes  equal  to  the  unnotched  fracture 
stress  as  the  flaw  size  is  reduced  to  zero.  Because  the  dcimage  zone  has  an 
effect  on  the  slot  crack  tip  stress  field  which  is  quantitatively  similar  to 
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TABLE  IV.  COMPARISON  OF  MEASURED  TENSILE  STRENGTH  OF  SEVERAL  LIMESTONES 
WITH  VALUES  CALCULATED  FROM  IHE  GRIFFITH  THEORY 
(Data  of  Perkins  and  Bartlett) 


Limestones 


Young's  Modulus 
GN/m2 


Tensile  Strength.  MN/m2 
Fracture  Energy  Measured,  Calculated*, 
R»  J/m2  MN/m2  MN/m2 


Austin 

11.0 

14.0 

1.81 

7.0 

Carthage 

32.4 

77.0 

4.68 

28.2 

Lueders 

19.3 

38.5 

4.11 

15.4 

Indiana 

19.3 

84.0 

3.58 

22.7 

*  Assuming  a  critical  flaw  size  of  1  mm. 
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FIGURE  15.  THE  GENERAL  RELATION  BETWEEN  FRACTURE  STRESS  AND  CRACK  LENGTH  IN 
MATERIALS  CAPABLE  OF  INELASTIC  ’RELAXATIONS  AT  THE  CRACK  TIP.  For 
large  crack  lengths  Equation  (A)  applies.  For  small  crack  lengths, 
as  at  aQ,  inelastic  deformation  becomes  predominant  and  the  dis¬ 
crepancy  between  the  predicted  fracture  stress,  at  (A),  and  the 
measured  fracture  stress,  at  (5) ,  can  be  quite  largeT 
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that  of  0  plastic  zone,  we  can  also  expect  a  very  similar  crack  length- 
fracture  stress  relation  in  rocks.  Thus,  if  a  rock  sample  contains  a  small 
flaw,  a  large  portion  of  the  rock  will  undergo  microcracking  before  the 
stress  on  the  rock  is  sufficient  to  extend  the  flaw.  In  such  a  case  a 
description  of  the  critical  condition  for  crack  extension  must  take  these 
inelastic  processes  into  account.  The  use  of  Equation  (4),  for  example,  to 
compute  R  from  the  unnotched  tensile  strength  and  an  assumed  flaw  size,  say, 

*  characteristic  of  the  microstructure  of  '.he  rock,  presupposes  a  fracture 
stress -flaw  size  relation  as  indicated  by  the  dashed  line  in  Figure  14.  The 
result  is  a  gross  underestimate  of  the  true  R. 

An  important  consequence  of  this  argument  is  the  expectation  that  R- 
values  of  rocks,  like  the  toughness  values  of  ductile  metals,  are  systematically 
underestimated  and  display  a  flaw-  and  specimen-size  dependence  when  these 
dimensions  are  relatively  small,  e.g.,  when 

D  >  C  E  R/ac2  ,  (5) 

where  D  is  the  limiting  specimen  dimension,  C  is  a  constant  of  order  0.5 
to  2  depending  on  the  specimen  geometry,  and  is  the  stress  at  which  micro¬ 
crack  generation  becomes  significant,  i.e.,  the  tensile  strength  of  an  un¬ 
notched  test  bar.  The  critical  specimen  dimensions  are  usually  the  specimen 
width  in  the  case  of  tensile  specimens  or  the  beam  height  in  the  case  of  a 
DCB  and  notch  bend  specimen.  For  values  of  E,  R,  and  O typical  for  limestone 
D  ~  10  cm  indicating  that  the  size  of  our  specimens  may  have  been  too  small  to 
measure  the  size -independent  value  of  R.  This  size  effect  is  distinct  from 
that  associated  with  Weibull-type  distributions  of  flaws  or  other  microstruc- 
tural  entities,  and  could  account  for  the  lower  R-value  reported  by  Friedman 
and  coworkers  for  limestone  (see  Table  II) . 
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It  should  be  apparent  that  the  microcrack -zone  of  damage  concept 
offers  new  opportunities  for  explaining  the  mechanical  behavior  of  brittle 
materials.  In  the  case  of  rocks  it  makes  it  possible  to  rationalize: 

(1)  The  discrepancy  between  R  and  surface  energy 

(2)  The  discrepancy  between  measured  and  calculated  unnotched 
strength  values 

(3)  The  lack  of  a  correlation  between  fracture  energy  and 
unnotched  strength  values  among  different  types  of  rocks 
and  ceramic  materials 

(4)  A  wide  range  of  microstructural  effects 

(5)  Size  effects  and 

(6)  The  apparent  insensitivity  of  R  to  notch  root  radius. 

Since  the  fracture  energies  of  some  ceramics  approach  those  of  rocks,  the 
microcrack-zone  of  damaga  approach  may  also  be  of  some  value  for  interpret¬ 
ing  and  modifying  the  strength  and  toughness  or  this  class  of  materials. 
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APPENDIX  A 


Computation  of  En»rgy  Dissipation  Rates 


At  the  basis  of  the  Griffith  energy  balance  or  thermodynamic  approach 

to  fracture  Is  the  calculation  of  the  stored  elastic  strain  energy  in  the  body 
given  by 


(A-D 

Here  c  denotes  elastic  strain  and  a  the  stress.  The  subscripts  follow  the 
Einstein  tensor  notation.  The  outer  integral,  r,  denotes  integration  over  the 
volume  of  the  body.  The  evaluation  of  U  is  greatly  simplified  if  there  exists 
a  linear  (Hcokeian)  relation  between  stress  and  strain  throughout  the  body. 

Further,  the  elastic  moduli  are  assumed  constant  throughout  and  that 
the  material  remains  compatible  (i.e.,  no  cracks  or  voids  develop). 

This  latter  assumption  is  apparently  violated  to  various  degrees  in  rocks 
subject  to  tensile  stresses.  As  a  consequence  of  forming  or  opening  many  small 
cracks,  the  material  becomes  more  compliant.  Equivalently,  the  bulk  elastic 
properties  appear  to  be  reduced .  However,  within  the  framework  of  these  assump¬ 
tions,  we  can  obtain  in  the  case  of  a  DCB  specimen,  a  simple  relation  between 
the  component  of  force  perpendicular  to  the  crack,  P,  and  displacement,  y  (between 
the  points  of  load  application) : 


where  E  is  Young's  modulus,  b  is  specimen  thickness,  and  0  is  the  compliance 

which  depends  only  on  the  specimen  dimensions  in  the  plane  perpendicular  to 

the  thickness  and  crack  length.  The  stored  elastic  strain  energy  also  simplifies 
to 

f\  -  I 


A-2 


u  ~  \  Py 


(A-3) 


and  substituting  from  Equation  (A-2) 


U 


(A-4) 


The  stored  elastic  energy  decreases  as  the  crack  extends  at  a  rate  given  by 


or 


(A-5) 


where  G  is  the  elastic  energy  release  rate.  If  the  stored  elastic  energy  is 
released  at  the  same  rate  as  it  is  dissipated  then  G  -  R  and  Equation  (A-5)  is 
the  means  of  calculating  R  provided  0  can  be  related  to  crack  length.  Kanninen^ 
has  derived  this  relation  and  obtains 


« 


-  -f-r  [ 2 \3a3  +  6A2a2  f S^-*c-C-°°h  ±  sln  \c  cos  As\  + 

X  h  L  \  Sinh  Xc  -  sin2  Xc  / 


+  6Xa  sinhlAc 


iSinh  Xc 


■f  sin2  Xc\ 
-  sin2  Xc/ 


+  3  (  Sinh  Xc  Cosh  Xc  -  sin  Xc  cos  Xc 
Sinh2  Xc  -  sin2  Xc 


(A-6) 


SSL 


1/4 


c  ■  L  -  a 


where  h  is  the  beam  height  (one-half  the  width  of  the  specimen)  and  L  is 
the  specimen  length. 
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Substituting  Equation  (A-6)  Into  (A-5)  give.  .  mcan.  of  determl.Jlng 

R  in  tanas  of  th,  assured  displacement  y  and  crack  length,  a.  The  displacement 
1.  measured  at  the  end  of  the  steel  pin  block.  (cf„  Plg„re  1).  Slnce  the  pin 
block,  hay,  negligible  compliance  compared  to  the  rock,  the  point  of  load  appli¬ 
cation  is  taken  as  the  Interface  between  the  end  of  the  specimen  and  the  pin 

blocks.  The  measured  displacement  1.  corrected  to  correspond  to  the  displacement 
at  this  point. 

A.  the  crack  extension  was  often  difficult  to  measure,  the  crack  length 
could  be  obtained  Indirectly  baaed  on  the  assumption  that  th.  coefficient  of 
friction  between  the  wedge  and  pins  remained  essentially  constant  during  the 
test.  We  have  P  -^Frf,ere  *  1.  the  compliance  relating  y  to  F,  the  force 
applied  to  the  wedge,  and  the  above  assumption  permit.  to  ..^constant . 

This  ratio  1.  obtained  at  the  beginning  of  the  test  from  the  slope  of  the  linear 
portion  of  the  F-y  curve  (cf„  Figure  2)  and  0  which  1.  calculated  from  Equation 
(A-6)  using  the  Initial  .lot  length.  For  other  value,  of  y  and  F,  the  compliance, 
0,  and  crack  length  can  he  obtained  from  the  following  combination  of  equations: 


w!  I 

*  i  v 

i  (t) 


0  ■  Aa 


(A-7) 


where  A  and  n  are  constant..  Equation  (A-7)  1,  empirical  but  ha.  been  shown  to 

provide  a  good  correlation  between  the  compliance  and  crack  length/9’  The 

constant,  A  and  „  are  obtained  from  a  least  square,  fit  of  Equation  (A-6)  to 
Equation  (A-7) . 
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Estimation  of  Total  Microcrack  Area 

Referring  to  Figure  B-l  we  can  estimate  the  eff  ct  of  a  regular 
array  of  microcracks  on  the  effective  Young's  modulus  of  a  material  by  assum¬ 
ing  that  the  cracks  are  noninteracting,  i.e.,  they  are  sufficiently  far  apart 
that  th-j  stress  field  around  any  crack  is  essentially  independent  of  the 
presence  of  neighboring  cracks.  The  compliance  of  a  rectangular  plate  of 
dimensions  ^  x  w  and  containing  one  crack  of  length  2a  is  given  approximately 
by 

*  '  (1  +  £T>  (»-» 

where  b  ie  the  thickness.  If  we  now  assemble  a  plate  of  dimensions  L  x  W 
from  tnese  singly  cracked  plates,  the  compliance  of  the  assembly  becomes 


4> 


L  t-i 
EbW  (1  + 


Pffcr2 
b  ' 


(B-2) 


In  Equation  (B-2),  p  is  the  number  of  microcracks  per  unit  volume,  given  by 
p  “  1/w-fcb  and  a  ■  ba,  the  area  per  microcrack.  Equation  (B-2)  is  strictly 
applicable  only  to  through-the -thickness  cracks  but  for  the  sake  of  simpli¬ 
fication,  we  apply  it  to  the  case  of  a  three-dimensional  array  of  microcracks 
each  completely  embedded.  Therefore,  by  comparison  with  the  compliance  of  an 
uncracked  plate  we  obtain 


(B-3) 


for  the  effective  Young's  modulus  of  a  solid  containing  imbedded  microcracks. 


For  a  DCB  specimen  under  load  it  can  be  recognized  that  if  the  load 
is  sufficiently  large  to  produce  microcracking  Egff  will  vary  throughout  the 


THE  MODEL  USED  TO  PREDICT  THE  EFFECT  OF  FLAW  DENSITY  ON 
COMPLIANCE.  The  plate  of  dimension  L  x  W  is  composed 
of  elements  of  dimension  i  x  w  each  containing  a  crack 
of  length  2a. 
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specimen  being  smallest  near  the  main  crack  tip  where  the  microcrack  density, 
and,  hence,  p,  is  the  greatest.  To  establish  a  reasonable  distribution  of 
microcracks  in  the  specimen  we  make  the  assumption  that  the  rock  has  a  con¬ 
stant  failure  stress,  cr^,  i.a.,  microcracking  begins  when  the  tensile  stress 

is  increased  to  C  but  can  become  no  larger  than  O’  .  This  assumption  allows 
c  c 

us  to  define  the  size  and  shape  of  the  damage  zone.  A  further  simplification 
is  obtained  by  considering  the  DCB  specimer  as  two  cantilever  beams  attached 
to  a  rigid  foundation.  Therefore,  as  each  beam  is  loaded  they  first  behave 
perfectly  elastic  until  0£  is  reached  at  the  point  of  maximum  stress  on  the 
tensile  side  of  the  beam.  At  any  stage  in  the  growth  of  the  damage  zone  we 

can  deduce,  from  application  of  inelastic  beam  theory,  that  the  effective 
Young's  modulus  at  any  point  within  the  damage  zone 


9Eh(|  - 


Jeff 


(B-4) 


9h(j  -  “)  +  2(C  -  y) 


where  S  =  h^b^/P 

x  =  coordinate  along  the  length  of  the  beam  measured 
from  the  foundation* 

y  =  coordinate  perpendicular  to  x  arid  measured  from 
the  tensile  fiber  of  the  beam* 

Q  =  the  value  of  y  at  the  damage  zone  boundary. 

Comparing  Equations  (B-4)  and  (B-3)  indicates  that 


b 


_  2(fi, -  Xl 


(B-5) 


*  The  extent  of  the  damage  zone  is  given  by  0  v  x  <  a  -  -  and  0  <  y  <  C 

3h(a-x)  h 
S  '  2  * 


where  C 
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If  we  note  that  the  term  pot  is  the  amount  of  microcrack  area  per  unit  volume, 
then  the  total  microcrack  area  in  the  damage  zone,  A,  is  obtained  by  integrat¬ 
ing  the  right  side  of  Equation  (B-5)  over  the  volume  of  the  damage  zone  giving 
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